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Abstract The f-sum rule is introduced and its applications to electronic and vi- 
brational modes are discussed. A related integral over the intra-band 
part of (j(tj) which is also valid for correlated electrons, becomes just 
the kinetic energy if the only hopping os between nearest neighbor sites. 
A summary is given of additional sum rule expressions for the optical 
conductivity and the dielectric function, including expressions for the 
first and second moment of the optical conductivity, and a relation be- 
tween the Coulomb energy and the energy loss function. It is shown 
from various examples, that the optical spectra of high Tc materials 
along the c-axis and in the ab-plane direction can be used to study the 
kinetic energy change due to the appearance of superconductivity. The 
results show, that the pairing mechanism is highly unconventional, and 
mostly associated with a lowering of kinetic energy parallel to the planes 
when pairs are formed. 

Keywords: Optical conductivity, spectral weight, sum rules, reflectivity, dielec- 
tric function, inelastic scattering, energy loss function, inelastic elec- 
tron scattering, Josephson plasmon, multi-layers, inter-layer tunneling, 
transverse optical plasmon, specific heat, pair correlation, kinetic en- 
ergy, correlation energy, internal energy. 



1. Macroscopic electromagnetic fields in matter 

1.1 Introduction 

The response of a system of electrons to an externally applied field 
is commonly indicated as the dielectric function, or alternatively as the 
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optical conductivity. The discussion in this chapter is devoted to induced 
currents and fields which are proportional to the external fields, the so- 
called linear response. The dielectric and the optical conductivity can be 
measured either using inelastic scattering of charged particles for which 
usually electrons are used, or by measuring the absorption of light, or 
the amplitude and/or phase of light reflected or transmitted by a sample. 
The two cases of fast particles and incident radiation involve different 
physics and will be discussed separately. 

1.2 Reflection and refraction of electromagnetic 



Optical spectroscopy measures the reflection and refraction of a beam 
of photons interacting with the solid. A rarely used alternative is the use 
of bolometric techniques to measure the absorption of photons directly. 
A variety of different experimental geometries can be used, depending on 
the type of sample under investigation, which can be a reflecting surface 
of a thick crystal, a free standing thin film, or a thin film supported by a 
substrate. Important factors influencing the type of analysis are also the 
orientation of the crystal or film surface, the angle of incidence of the 
ray of photons, and the polarization of the light. In most cases only the 
amplitude of the refiected or refracted light is measured, but sometimes 
the phase is measured, or the phase difference between two incident rays 
with different polarization as in ellipsometry. The task of relating the 
intensity and/or phase of the refiected or refracted light to the dielectric 
tensor inside the material boils down to solving the Maxwell equations 
at the vacuum/sample, sample/substrate, etc. interfaces. An example is 
the ratio of the reflection coefficients {Rp/ Rg jand phase differences (?7p — 
rjs) of light rays with p and s-polarization reflected on a crystal/vacuum 
interface at an angle of incidence 9. These quantities which are measured 
directly using ellipsometry 



The real and imaginary part of the dielectric constant can be calculated 
from such a measurement with the aid of Eq. 1.1. In contrast to a beam 
of charged particles, the electric fleld of a plane electromagnetic wave 
is transverse to the photon momentum. The dielectric tensor elements 
which can be measured in an optical experiment are therefore transverse 
to the direction of propagation of the electromagnetic wave. In a typical 
optical experiment the photon energy is below 6 eV. In vacuum the 
photon wave number used in optical experiments is therefore O.OOOSA"^, 
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or smaller, which is at least three orders of magnitude below the Fermi 
momentum of electrons in a solid. For this reason it is usually said that 

optical spectroscopy measures the transverse dielectric constant or the 
optical conductivity at zero momentum. The optical conductivity tensor 
expresses the current response to an electrical field 

J(f, t)= j d'r' J dt'a{f, r', t - t'Mr', t') (1.2) 

Prom the Maxwell equations it can be shown that for polarization trans- 
verse to the propagation of an electromagnetic wave dE/dt = dD/dt + 
47rj. If the sample has translational invariance, the optical conductivity 
tensor has a diagonal representation in fe-space. Due to the fact that the 
translational symmetries of a crystalline solid are restricted to a discrete 
space group, k is limited to the first Brillouin zone. Consequently, as 
shown by Hanke and Sham[l] , the /c-space representation of the dielectric 
tensor becomes a matrix in reciprocal space 

a{G,u)g^^, = j d^rj rfV j dte'^^+^>^e-'^^+^'>^e'''*a{ry,t) (1.3) 

The dependence of a{q,u))g ^, on the reciprocal lattice vectors G, G' 
reflects, that the local fields can have strong variations in direction 
and magnitude on the length scale of a unit cell. Yet due to the long 
wavelength of the external light rays the Fresnel equations involve only 
G = G' = 0. Usually in texts on optical properties the only optical 
tensor elements considered have G = G' = 0, and in this chapter we will 
do the same. 

Inside a solid the wavelength of the electromagnetic rays can be much 
shorter than that of a ray with the same frequency travelling in vacuum. 
Although in this chapter we will not encounter experiments where the 
finite momentum of the photon plays an important role, we should keep 
in mind that in principle the photon momentum is non-zero and can have 
a non-trivial effect on the optical spectra. In particular it may corrupt 
Kramers-Kronig relations, which is just one out of several reasons why 
spectroscopic ellipsometry should be favored. 

1.3 Inelastic scattering of charged particles 

When a fast charged particle, moving at a velocity Ue, interacts weakly 
with a solid, it may recoil inelastically by transferring part of its mo- 
mentum, hq and its energy, fiu) to the solid. The fast electron behaves 
like a test charge of frequency lo = q- Vg, which corresponds to a dielec- 
tric displacement field, D{f,t) = eq~^ exp {iq ■ r — iut) . The dielectric 
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displacement of the external charges may be characterized by a density 
fluctuation, which has no field component transverse to the wave. D(f, t) 
is therefore a purely longitudinal field. In a solid mixing of transverse 
and longitudinal modes occurs whenever fields propagate in a direction 
which is not a high symmetry direction of the crystal. However, in the 
long wavelength limit the dielectric properties can be described by only 
three tensor elements which correspond to the three optical axes of the 
crystal. Since along these directions no mixing between longitudinal and 
transverse response occurs, we will consider the situation in this chapter 
where the fields and their propagation vector point along the optical 
axis. Inside a material the dielectric displacement is screened by the 
response of the matter particles, resulting in the screened field E{r, t) 
inside the solid[l]. 

E{f,t) = J d^f' J dt'e-\f,f',t-t')D(f',t') (1.4) 

For the same reasons as for the optical conductivity the A;-space repre- 
sentation of the dielectric tensor becomes a tensor in reciprocal space 

t-\q+G',q+G,u) = Jd^rJ d^f^ J dte'^'i+^'>''e-'^'i+^>'"' e'^^'e-^r^ ,t) 

(1.5) 

where G and G' denote reciprocal lattice vectors. The relation between 
the dielectric displacement and the electric field is 

E{q + G', io)=J2 + G', q + G, u)D{q + G, co) (1.6) 

G 

The macroscopic dielectric constant, which measures the macroscopic 
response to a macroscopic perturbation, i.e. for vanishingly small q, is 
given by[l] 

e{Lo) = lim — T -i ^ — r (1.7) 
q-.Oe-^{q,q,uj) 

where it is important, that in this expression first the matrix e{q + 
G', q + G^uj) has to be inverted in reciprocal space, and in the next step 
the (G = 0, G' = 0) matrix element is taken of the inverted matrix 
[1]. Energy loss spectroscopy using charged particles can be used to 
measure the dielectric response as a function of both frequency and 
momentum. This technique provides the longitudinal dielectric function, 
i.e. the response to a dielectric displacement field which is parallel to 
the transferred momentum q. The probability per unit time that a 
fast electron transfers momentum q and energy hu to the electrons was 
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derived by Nozieres and Pines [2, 3] for a fully translational invariant 
'jellium' of interacting electrons 

P(g,c.) = ^Im|— i-j (1.8) 

where e is the elementary charge. 

1.4 Relation between (t{(jj) and e{(jj) 

We close this introduction by remarking, that for electromagnetic 
fields propagating at a long wavelength the two responses, longitudi- 
nal and transverse, although different at any nonzero wave vector, are 
very closely related. We will take advantage of this fact when later in 
this chapter we extract the energy loss function for g « from optical 
data. According to Maxwell's equations for g ^ the uniform cur- 
rent density is just the time derivative of the uniform dipole field, hence 
iirj = iuj{E — D). Consequently for g ^ the conductivity and the 
dielectric function are related in the following way 

A/jri. 

e(0,a;) = l + (T(0,a;) (1.9) 

Throughout this chapter we will use this identity repeatedly. 

2. Interaction of light with matter 
2.1 The optical conductivity 

Let us now turn to the discussion of the microscopic properties of 
the optical conductivity function. The full Hamiltonian describing the 
electrons and their interactions is 

H = Y.^Aa<^ka + Y.UGP-G + \Y.ymp-k (1-10) 

ka G k 

k = (1.11) 

In this expression the symbol creates a plane wave of momentum Up 

and spin-quantum number cr, Uq represents the potential landscape due 
to the crystal environment. The third term is a model electron-electron 
interaction Hamiltonian, representing all electron-phonon mediated and 
Coulomb interactions, where is the fe'th Fourier component of the 
density operator. In addition to the direct Coulomb interaction, various 
other contributions may be relevant, such as direct exchange terms. As a 
result the spin- and momentum dependence of the total interaction can 
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have a more complex form than the above model Hamiltonian. Relevant 
for the subsequent discussion is only, that the interaction term commutes 
with the current operator. The quantum mechanical expression for the 
current operator is 

3? = E ^^i-g/2,a W2,a (1-12) 
p,cr 

The current and density operators are symmetric in k-space, satisfying 
pI. = p-k and = j-k- In coordinate space the representations of the 
density and the current are 

^(^1 = ^E«^''% (1-13) 

= ^Y.^''% (1-14) 

It is easy to verify, that together n(r) and j{r) satisfy the continuity 
equation iefi'^[n{r), H] + V ■ j{r) = 0. 

Let us now consider a many-body system with cigcn-statcs \m) and 
corresponding energies Em- For such a system the microscopic expres- 
sion for the optical conductivity has been explained by A.J. Millis in 
Chapter 6. The result for finite q was derived in 'The theory of quan- 
tum liquids' part I, by Nozieres and Pines (equation 4.163). For brevity 
of notation we represent the matrix elements of the current operators as 

i2:^^(n|ia,Jm) (1.15) 

With the help of these matrix elements, and with the definition Uumn = 
Em — En the expression for the optical conductivity is 

• 2 AT ■ J3(n-E„) r Anmjmn jnm jmn 

^ , .\ — I Ja,qJa,-q Ja,-qJa,q 

'^a,a{Q,'^ — — TP r > — — ; 

mVu ^-^^ Vu ItU-Umn + iV U + Umn + 

(1.16) 

Here N is the number of electrons, V the volume, m the electron mass, Qe 
the elementary charge, O is the thermodynamic potential, P = l/ksT 
and r/ is an infinitesimally small positive number. In principle in the 
calculation of Eq. 1.16 terms may occur under the summation for which 
^mn = 0. As (jjjnn occurs in the denominator of this expression, these 
zeros should be cancelled exactly by zeros of the current matrix elements, 
which poses a special mathematical challenge. 

In Eq. 1.16 (j{oj) is represented by two separate terms, a 5-function 
for a; = and a summation over excited many-body eigen-states. The 
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(5-function is a diamagnetic contribution of all electrons in the system, 
the presence of which is a consequence of the gauge invariant treatment 
of the optical conductivity, as explained by Millis in chapter 6. The 
presence of this term is at first glance rather confusing, since left by 
itself this 5-function would imply that all materials (including diamond) 
are ideal conductors! However, the second term has, besides a series 
of poles corresponding to the optical transitions, also a pole for = 0, 
corresponding to a negative 5- function of Recj(ix'). It turns out, that for 
all materials except ideal conductors this (5-function compensates exactly 
the first (diamagnetic) term of Eq. 1.16. This exact compensation is a 
consequence of the relation [6] 

■nm -mn at 2 

For every n: V ^-^£2^ = tl^ (1.17) 



Experimentally truly 'ideal' conductivity is only seen in superconduc- 
tors. In ordinary conducting materials the diamagnetic term broadens 
to a Lorentzian peak due to elastic and/or inelastic scattering. The 
width of this peak is the inverse life-time of the charge carriers. Of- 
ten in the theoretical literature the broadening is not important, and 
the Drude peak is counted to the Dirac-function in the origin. The in- 
frared properties of superconductors are characterized by the presence 
of both a purely reactive diamagnetic response, and a regular dissipa- 
tive conductivity [5]. The sum of these contributions counts the partial 
intra-band spectral weight which we discuss in section 2.4. With the 
help of Eq.1.17, the diamagnetic term of Eq. 1.16 can now be absorbed 
in the summation on the right-hand side 

a,.(g,cc;) = i V I Jc.,<iJa,-^ ^ Ja,-<iJa,g I 

(1.18) 

As explained in section 1.2, usually in optical experiments one assumes 
q^Qin the expressions for a{oj). It is useful at this stage to introduce 
the generalized plasma frequencies = '&-Kef^^^~^^^\j'^\^ui:^^V~^ , 

with the help of which we obtain the following compact expression for 
the optical conductivity tensor 

^aa(a;) = |^ E , (1-19) 

Although formally the parameter is understood to be an infinitesi- 
mally small positive number, a natural modification of Eq.1.19 consists 
of limiting the summation to a set of oscillators representing the main 
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optical transitions and inserting a finite value for jmn, which in this 
case represents the inverse lifetime of the corresponding excited state 
(e.g. calculated using Fermi's Golden Rule). With this modification 
Eq.1.19 is one of the most commonly used phenomenological representa- 
tions of the optical conductivity, generally known as the Drude-Lorentz 
expression. 

2.2 The f-sum rule 

The expressions Eqs. 1.16, 1.18, and 1.19 satisfy a famous sum rule. 
This is obtained by first showing with the help of Eq. 1.17, that for each 
n 

E - ^e/^(---") (1.20) 

m,m^n 

Second, as a result of Cauchy's theorem in Eq.1.19 the integral over all 
(positive and negative) frequencies of fRea{Lo) equals J^mn^mn/'^- To 
complete the derivation of the f-sum rule we also use that J2n e^(^~^") = 
1, which follows from the definition of the thermodynamic potential. 
Then 



J — ( 



Rea(u)du; = (1.21) 
, mV 



is the f-sum rule, or Thomas-Reich-Kuhn rule. It is a cornerstone for 
optical studies of materials, since it relates the integrated optical conduc- 
tivity directly to the density of charged objects, and the absolute value 
of their charge and mass. It reflects the fundamental property that also 
in strongly correlated matter the number of electrons is conserved. Note 
that the right-hand side of the f-sum rule is independent of the value of 
fi. Also the f-sum rule applies to bosons and fermions alike. Because 
Rea"(a;)=Rc(T(— a;) the sum rule is often presented as an integral of the 
conductivity over positive frequencies only. Superconductors present a 
special case, since Rccr(aj) now has a 5 function at = 0: Only half of 
the spectral weight of this (5-function should be counted to the positive 
frequency side of the spectrum. 

2.3 Spectral weight of electrons and optical 
phonons 

The optical conductivity has contributions from both electrons and 
nuclei because each of these particles carries electrical charge. The inte- 
gral over the optical conductivity can then be extended to the summation 
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over all species of particles in the solid with mass rrij, and charge qj 



Because the mass of an electron is several orders of magnitude lower than 
the mass of a proton, in many cases the contribution of the nuclei to the 
f-sum rule is ignored in calculations of the integrated spectral weight of 
metals. However, important exceptions exist where the phonon contri- 
bution cannot be neglected, notably in the c-axis response of cuprate 
high Tc superconductors. Although Eq.1.22 is completely general, in 
practice it cannot be applied to experimental spectra directly. This is 
due to the fact that the contributions of all electrons and nuclei can only 
be obtained if the conductivity can be measured sufficiently accurately 
up to infinite frequencies. In practice one always uses a finite cut-off. 
Let us consider the example of an ionic insulator: If the integral is car- 
ried out for frequencies including all the vibrational modes, but does not 
include any of the inter-band transitions, then the degrees of freedom 
describing the motion of electrons relative to the ions is not counted. As 
a result the large number of electrons and nuclei which typically form 
the ions are not counted as separate entities. Effectively the ions behave 
as the only (composite) particles in such a case, and the right-hand side 
of Eq. 1.22 contains a summation over the ions in the solid. Application 
of Eq.1.22 provides the so-called transverse effective charge, which for 
ionic insulators with a large insulating gap corresponds rather closely 
to the actual charge of the ions. In the top panel of Fig. 1.1 this is 
illustrated with the infrared spectrum of MgO. Indeed the transverse 
effective charge obtained from the sum rule is 1.99, in good agreement 
with the formal charges of the Mg^"^ and ions. 

Because the mass of the ions is much higher than the free electron 
mass, the corresponding spectral weight integrated over the vibrational 
part of the spectrum is rather small. In a metal, even if optical phonons 
are present, usually the spectral weight at low frequencies is completely 
dominated by the electronic contributions due to the fact that the free 
electron mass is much smaller than the nuclear mass. A widely spread 
misconception is, that the screening of optical phonons in metals leads 
to a smaller oscillator strength than in ionic insulators. The opposite 
is true: Due to resonant coupling between vibrational modes and elec- 
tronic oscillators, the optical phonons in an intermetallic compound of- 
ten have much more spectral weight than optical phonons in insula- 
tors. This 'charged phonon' effect was formulated in an elegant way 
in 1977 by Rice, Lipari and Strassler[9]. They demonstrated, that un- 
der resonant conditions, due to electron-phonon coupling, vibrational 
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Figure 1.1. Optical conductivity of MgO (top panel) and FeSi at r = 4 K (bottom 
panel). In the insets the function Z{u)) = 8/i(47rnie^)~°'^ crph(ai')dw'} ' is dis- 
played. For FeSi the electronic background (dotted curve of the lower panel) was 
subtracted. For uj > 600cm~^ •^('^) corresponds to transverse effective charge. Data 
from Ref. [7, 8]. 

modes borrow oscillator strength from electronic modes, which boosts 
the intensity of the vibrational modes in the optical conductivity spec- 
tra. This effect is now known to be common in many materials, for 
example in TCNQ-salts,blue bronze, IV- VI narrow-gap semiconductors, 
FeSi and related compounds, and the beta-phase of sodium vanadate 
[10, 11, 12, 13, 15, 8, 7, 14, 16]. In Figs. 1.1 and 1.2 the charged phonon 
effect is illustrated using the examples of FeSi, MnSi, CoSi and RuSi 
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[8, 7, 14], showing that the transverse charge is between 4 and 5. These 
compounds are not ionic insulators, because the TM and Si atoms have 
practically the same clcctro-ncgativities and electron affinities. Instead 
the large transverse charge of these compounds arises from the charged 
phonon effect predicted by Rice. The strong temperature dependence of 
the transverse charge of FcSi correlates with the gradual disappearance 
of the semiconductor gap as the temperature is raised from 4 to 300 K. 
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Figure 1.2. Temperature dependence of the transverse effective charge of Co-Si, Fe- 
Si, Ru-Si, and Mn-Si pairs, calculated from the oscillator strength of the optical 
phonons. Data from Refs. [7], [8] and [14]. 

2.4 Partial spectral weight of the intra-band 
transitions 

Often there is a special interest in the spectral properties of the charge 
carriers. The electrons are subject to the periodic potential of the nu- 
clei, resulting in an energy-momentum dispersion which differs from free 
electrons. Often one takes this dispersion relation as the starting point 
for models of interacting electrons. The Coulomb interaction and other 
(e.g. phonon mediated) interactions present the real theoretical chal- 
lenge. The total Hamiltonian describing the electrons and their interac- 
tions is then 

k,cr kqp ac' 
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The current operator is in this case 

^9 = I E - Cp_g/2,<,Cp+q/2,a (1-24) 

where the = h~^dek/dk is the group velocity. The density operator 
commutes with the interaction part of the Hamiltonian. This has an 
interesting and very useful consequence, namely that a partial sum rule 
similar to the f-sum rule exists, which can be used to probe experimen- 
tally the kinetic energy term of the Hamiltonian. This partial sum rule 
for integration of the intra-band conductivity, crj^ai'^)' yields[17] 

/oo 
-00 

1 

nik 

Apparently the total spectral weight contained in the inter-band transi- 
tions, (7^(0;), is exactly 

In the limit where the interaction Vfc = 0, the occupation function in 
the above summation is a step-function at the Fermi momentum. In this 
case the summation over k becomes an integration over the Fermi volume 
with Uk set equal to 1. After applying Gauss's theorem we immediately 
obtain the well-known Fermi surface integral formula 

^p,a = 9 Isp iva{a)daa (1.28) 

where g is the spin degeneracy factor. In the literature two limiting 
cases are most frequently considered: (i) the free electron approxima- 
tion, where = m,g is the free electron mass independent of the mo- 
mentum of the electron, and (ii) the nearest neighbor tight-binding limit. 
In the latter case the dispersion is = —2tx cos kxttx — 2ty cos kyUy — 
2tzCOskzaz, with the effect that l/rrik = —2thr'^a^cos{kaaa), and 

^ Keaac.{u:)dw = - Ek,a rikek = {-Hkin) (1-29) 

where the integration should be carried out over all transitions within 
this band, including the 5-function at w = in the superconducting 
state. The upper limit of the integration is formally represented by 
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the upper limit Um- In practice the cutoff cannot always be sharply 
defined, because usually there is some overlap between the region of 
transitions within the partially filled band and the transitions between 
different bands. Hence in the nearest neighbor tight-binding limit the 
f-sum provides the kinetic energy contribution, which depends both on 
the number of particles and the hopping parameter t[18]. This relation 
was used by Baeriswyl et al. to show, using exact results for one di- 
mension, that the oscillator strength of optical absorption is strongly 
suppressed if the on-site electron-electron interactions (expressed by the 
Hubbard parameter U) are incrcascd[19]. The same equation can also 
be applied to superconductors, examples will be discussed later in this 
chapter. In the case of a superconductor it is important to realize, that 
the integration on the left-hand side of Eq. 1.29 should also include the 
condensate (5- function at = 0. As the optical conductivity can only be 
measured for u; > 0, the spectral weight in the 5-function has to derived 
from a measurement of the imaginary part of C7(a;), taking advantage of 
the fact that the real and imaginary part of a 5-function conductivity 
are of the form 

^singular/, ,\ ^^^p,s 

The plasma frequency of the condensate, ujp^s, is inversely proportional 
to the London penetration depth, A(T) = c/ujp^s{T) with c the velocity 
of light. In the literature [20, 21] the (^-function, conductivity integral for 
a; > 0, and the kinetic energy are sometimes rearranged in the form 

4" = "^i-Hkin) - /or Rea(a;)da; (1-30) 

Whenever the kinetic energy term on the right-hand side changes its 
value, this expression suggests a 'violation' of the f-sum rule, since the 
spectral weight in the (^-function now no longer compensates the change 
of spectral weight in the conductivity integration on the right-hand side. 
Of course there is no real violation, but part of the optical spectral weight 
is being swapped between the intra-band transitions and the inter-band 
transitions. Later in this chapter we will use the relation between kinetic 
energy as expressed in the original incarnation due to Maldague [18] (Eq. 
1.29) to determine in detail the temperature dependence of the a6-plane 
kinetic energy of some of the high Tg superconductors. 

It is easy to see, that for a small filling fraction of the band Eq. 1.29 is 
the same as the Galilean invariant result: The occupied electron states 
are now all located just above the bottom of the band, with an energy 
—t. Hence in leading orders of the filling fraction —{tpg\Ht\'il)g) = Nt. 
Identifying a??i~H~^ as the effective mass m* we recognize the familiar 
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f-sum rule, Eq.1.22, with the free electron mass replaced by the effective 
mass. 

As the total spectral weight (intra-band plus inter-band) should sat- 
isfy the f-sum rule, the intra-band spectral weight is bounded from above, 
i.e. < Yl]^nk/mk < n/m. Near the top of the band the dynamical 
mass has the peculiar property that it is negative, < 0, which in the 
present context adds a negative contribution to the intra-band spectral 
weight. On the other hand, the fact that Re(T(a;) has to be larger than 
zero, implies that the equilibrium momentum distribution function nfc 
is subject to certain bounds: If for example would preferentially oc- 
cupy states near the top of the band, leaving the states at the bottom 
empty, the intra-band spectral weight would acquire an unphysical neg- 
ative value. Apparently such momentum distribution functions cannot 
result from the interactions of Eq.1.23, regardless of the strength and 
k-dependence of those interactions. 

2.5 Additional sum rules for cr(u>) and l/e(u>) 

Several other sum rule type expressions exist for the optical conduc- 
tivity and for the dielectric constant. Here we give a summary. In the 
presence of a magnetic field an optical analogue of the Hall effect exists. 
The behavior is similar to the DC-limit, resulting in an off-diagonal com- 
ponent of the optical conductivity axy{u)) = —ayx{uo). where the z-axes 
is parallel to the magnetic field. The optical Hall angle is defined as 

tH{u) = ^ (1.31) 

The optical (axx) and Hall conductivities((Ta;y) can be measured directly 
in optical transmission experiments [22, 23]. Drew and Coleman have 
shown [24] that this response function obeys the sum rule 

2 

- / tH{uj)dw = ujH (1.32) 
vr Jo 

where the Hall frequency loh is unaffected by interactions, and in the 
Galilean invariant case corresponds to the bare cyclotron frequency, 

UJH = eB/m. 

A first moment sum rule of the optical conductivity is easily obtained 
for T = 0, by direct integration of Eq. 1.18, providing 

In free space there is no scattering potential nor a periodic potential 
causing Umklapp scattering. Hence for electrons moving in free space 
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the right-hand side of Eq. 1.33 is exactly zero. This comes as no sur- 
prise: The integral on the left-hand side is also zero, since the optical 
conductivity of such a system has only a (5-function at (j = due to 
Galilean invariance. However, in the presence of Umklapp scattering 
the eigen-states of the electrons with energy-momentum dispersion 
are no longer the free electron states in the summation of Eq. 1.33. The 
true eigen-states are superpositions of plane waves. Vice versa the free 
electron states generated by the c|. operators of the above expression 
can be written as a superposition of the eigen-states of the periodic po- 
tential: c\_^Q ^ = Z^m "^(3(^)^1 m cr' where the latter operator generates 
the m'th eigen-state with momentum k in the first Brillouin zone. For 
brevity we introduce the notation Aq = |a^(A;)p, and = o}^ j cr(^k,j,a- 
Expressed in terms of these band occupation number operators Eq. 1.33 
is 



/■oo ZTTP n 

hm / a;a«,«(g,u;)du; = ^- + G^f E ^g^gR(1 " O) 

^ k,G j,m,a 

(1.34) 

The summation on the right-hand side strongly suggests an intimate 
relationship between the optical conductivity and the kinetic energy of 
the electrons. However, due to the fact that the expression on the right- 
hand side is rather difficult to calculate, the first moment of a"(a;) is of 
little practical importance. It's main purpose in the present context is 
to demonstrate the trend that an increase of the kinetic energy is ac- 
companied by an increase of the first moment of the optical conductivity 
spectrum. This is consistent with the notion, that an increase of kinetic 
energy is accompanied by a blue-shift of the spectral weight. 

For the energy-loss function a separate series of sum rule type equa- 
tions can be derived[25, 26, 27] 



f 



Im— -dw = — (1.35) 

, e{Lij) mV 



which is similar to the f-sum rule for the optical conductivity, Eq. 1.22. 

As a result of the fact that the real and imaginary part of the energy 
loss function are connected via Kramers-Kronig relations, the following 
relation exists ^ 

Im — T^d^ = TT (1.36) 
-tx> u>e{u) 

This expression can in principle be used to calibrate the absolute inten- 
sity of an energy-loss spectrum, or to check the experimental equipment, 
since the right-hand side does not depend on any parameter of the ma- 
terial of which the spectrum is taken. We can use the relation between 
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e(a;) and a{io), Eq. 1.9, to express Eq. 1.36 as a function of a(Lo). Using 
Cauchy's theorem, it is quite easy to prove from Eq.1.36, that 

/ Re——-—dw = — (1.37) 

Often the intra-band optical conductivity is analyzed in terms of a fre- 
quency dependent scattering rate 1/t{lo) = {ne'^ /m)Re{a{uj)^^}, which 
follows directly from the experimental real and imaginary part of the 
optical conductivity. Taking the limit A ^ in Eq.1.37, we observe that 

1 , TT ne^ 



, du! = lim — = oo (1.38) 

/o ''"(i^) A^o 2A m 

Hence ultraviolet divergency appears to be a burden of integral formulas 
of the frequency dependent scattering rate[28, 29, 30] which is hard to 
avoid. 

In section 3.1 we will encounter a relation between the loss-function 
and the Coulomb energy stored in the electron fluid [26] 

This expression is limited to the ground state at T = 0, as was also the 
case for Eq. 1.33. The integrands on the left-hand side of Eq. 1.33 and 
Eq. 1.39 are odd functions of frequency. In contrast the f-sum rule, and 
the other expressions given in this subsection all involve integrals over 
an even function of frequency, which is the reason why the latter can be 
represented as integrals over all (positive and negative) frequencies. The 
fact that h occurs on the right-hand side of Eqs. 1.33 and 1.39 implies 
that these expressions are of a fundamental quantum mechanical nature, 
with no equivalent in classical physics. 

Recently Turlakov and Leggett derived an expression for the third 
moment of the energy loss function, which in the limit of —> is a 
function of the Umklapp potential of Eq.1.11 

/oo —UJ'^ Att"^ I \ 

-J^7^)'^ = ^{-^°''-"-^--) 

The fact, that the right-hand side of Eq. 1.40 is finite implies, that for 

a; ^ oo the loss function of any substance must decay more rapidly 
than Im{— €(a;)^^} oc w"^, and that the optical conductivity decays 
faster than Re{cr(ci;)} oc u!~^. This expression is potentially interesting 
for the measurement of changes in Umklapp potential, provided that 
experimental data can be collected up to sufficiently high photon energy, 
so that the left-hand side of the expression reaches its high frequency 
limit. 
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3. The internal energy of superconductors 

A necessary condition for the existence of superconductivity is, that 
the free energy of the superconducting state is lower than that of the 
non-superconducting state. At sufficiently high temperature important 
contributions to the free energy are due to the entropy. These contribu- 
tions depend strongly on the nature of the low energy excitations, first 
and foremost of all their nature be it fermionic, bosonic or of a more 
complex character due to electron correlation effects. At T = the free 
energy and internal energy are equal, and are given by the quantum 
expectation value of the Hamiltonian, which can be separated into an 
interaction energy and a kinetic energy. 

3.1 Interaction energy in BCS theory 



s-wave symmetry d-wave symmetry 
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Figure 1.3. The k-space representation of the superconductivity induced change of 

pair correlation function for the s-wave (left panel) and d-wave symmetry (right 
panel). Parameters: A/W = 0.2, lod/W = 0.2. Doping level x = 0.25 

We consider a system of electrons interacting via the interaction Hamil- 
tonian given in Eq. 1.23. In the ground state of the system, the inter- 
action energy, including the correlation energy beyond the Hartree-Fock 
approximation, is just the quantum expectation value of the second (in- 
teraction) term of 1.23. Here we are only interested in the difference in 
interaction energy between the normal and superconducting state. 

E'corr - Korr = T.^k {{PkP-k)s ' {PkP-k)n) = E ^^Sgk (1-41) 

k k 

In BCS theory the only terms of the interaction Hamiltonian which 
contribute to the pairing are the so-called reduced terms, i.e. those terms 
in the summation of Eq. 1.23 for which the center of mass momentum 
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p + q = 0. The quantum mechanical expectation value of the correlation 
function is 

Sgk = "^{lup+kl^ - Op+k)(.Op - \up\^) +'^Up+kVp+kulv* (1.42) 
p p 

The first term on the right-hand represents the change in exchange cor- 
relations, whereas the second term represents the particle-hole mixing 
which is characteristic for the BCS state. A quantity of special interest is 
the real space correlation function 5g{r, r') = {n{r)n{r'))s — {n{r)n{r'))n- 
The Fourier transform of this correlation function is directly related to 
5gk appearing in the expression of the interaction energy, Eq.1.41 

^9k = ^J d'r J d^'e'^^'-^'^Sgiry) (1.43) 

We see, that if the correlation function 6g{r, r') could be measured some- 
how, and the interaction is known, than the interaction energy would 
follow directly from our knowledge of Sg{r,r'): 

El„ - El',,, = Jd'rJ d\'V{r - r')5g{r, r') (1.44) 

In a conventional superconductor the quasi-particles of the normal state 
are also the fermions which become paired in the superconducting state. 
(Note, that now we are using the concept of Landau Fermi-liquid quasi- 
particles for the normal state. Later in this manuscript we will ex- 
plore some consequences of not having a Fermi liquid in the normal 
state, where the quasi-particle concept will be abandoned.) Although 
the quasi-particle eigcn-statcs of a conventional Fermi liquid have an 
amount of electron character different from zero, their effective masses, 
velocities and scattering rates are renormalized. The conventional point 
of view is, that pairing (enhancement of pair correlations) reduces the 
interaction energy of the electrons, by virtue of the fact that in the super- 
conducting state the pair correlation function g{r,r') = {'^\h{r)h{r')\'^) 
increases at distances shorter than the superconducting coherence length 
^0- If the interaction energy V{r — r') is attractive for those distances, 
the interaction energy, Eq. 1.44, decreases in the superconducting state, 
and V{r — r') represents a (or the) pairing mechanism. In Fig. 1.3 we 
show calculations of 5gk assuming a bandstructure of the form 

W 

Cfe = [cos kxa + cos kya] — fi (1-45) 
while adopting an order parameter of the form 



Ajfc = Aoe(|ejt - nl -ud) 



(1.46) 
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Figure 1.4- The coordinate space representation of the superconductivity induced 
change of pair correlation function for the s-wave (left panel) and d-wave symmetry 
(right panel). Parameters: A/W = 0.2, ujd/W = 0.2. Doping level: x = 0.25 

for S-wave symmetry, and 

Afe = Aq [cos k^a — cos kyo] Q{\ek — /x| — cjd) (1-47) 

for d-wave symmetry. The parameters used were A/W = 0.2, lod/W = 
0.2, and Ep/W = 0.43 corresponding to x=0.25 hole doping counted 
from half filling of the band. The chemical potential in the supercon- 
ducting state was calculated selfconsistently in order to keep the hole 
doping at the fixed value of x=0.25 [31, 32, 33, 34]. Prom Fig.1.3 we 
conclude that s-wavc pairing symmetry requires a negative regard- 
less of the value of k, whereas the d-wave symmetry can be stabilized 
either assuming > for fc in the (vr, vr) region, or Vfe < for k near 
the origin. Both types of symmetry are suppressed by having > at 
small momentum, such as the Coulomb interaction. 

In Fig. 1.4 we display the correlation function in coordinate space 
representation. This graph demonstrates, that d-wave pairing is stabi- 
lized by a nearest-neighbor attractive interaction potential. An on-site 
repulsion has no influence on the pairing energy, since the pair correla- 
tion function has zero amplitude for r — r' = 0. On the other hand, for 
s-wave pairing the 'best' interaction is an on-site attractive potential, 
since the s-wave 6g{r,r') reaches it's maximum value at r — r' = 0. 

3.2 Experimental measurements of the Coulomb 
interaction energy 

In a series of papers Leggett has discussed the change of Coulomb 
correlation energy for a system which becomes superconducting [35], and 
has argued, that this energy would actually decrease in the supercon- 
ducting state. Experimentally the changes of Coulomb energy can be 
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measured directly in the sector of /c-space of vanishing k. The best, 
and most stable, experimental technique is to measure the dielectric 
function using spectroscopic ellipsometry, and to follow the changes as 
a function of temperature carefully as a function of temperature. Be- 
cause the cuprates are strongly anisotropic materials, it is crucial to 
measure both the in-plane and out-of-plane pseudo-dielectric functions, 
from which the full dielectric tensor elements along the optical axes of 
the crystal then have to be calculated. We followed this procedure for 
a number of different high Tc cuprates, indicating that the Coulomb 
energy in the superconducting state increases for k=0. However, for 
k ^ this need no longer be the case. Summarizing the situation [36]: 
the Coulomb interaction energy increases in the superconducting state 
for small k. This implies, that the lowering of internal energy in the 
superconducting state must be caused either by other sectors of k-space 
(in particular at around the (vr, tt) point, see Fig. 1.3!), or by a lowering 
of the kinetic energy in the superconducting state. The latter is only 
possible in a non Fermi liquid scenario of the normal state. 

3.3 Kinetic energy in BCS theory 



1.0 ^ 1.0 r- 




kp k K k 



Figure 1.5. Occupation function as a function of momentum in tlie normal (dash) 
and the superconducting (soUd) state for Fermi Uquid (left panel) and an example of 
a broad distribution function, not corresponding to a Fermi liquid (right panel). 

In BCS theory the lowering of the pair-interaction energy is partly 

compensated by a change of kinetic energy of opposite sign. This can 
be understood qualitatively in the following way: The correlated motion 
in pairs causes a localization of the relative coordinates of electrons, 
thereby increasing the relative momentum and the kinetic energy of the 
electrons. Another way to sec this, is that in the superconducting state 
the step of at the Fermi momentum is smoothed, as indicated in the 
left panel of Fig. 1.5, causing E^in to become larger[37]. 
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A pedagogical example where the kinetic energy of a pair is higher 
in the superconducting state, is provided by the negative U Hubbard 
model[38]: Without interactions, the kinetic energy is provided by the 
expression 

Ekin = -t (*|c.W + i^.c.|*) (1.48) 

<i,j>,(T 

Let us consider a 2D square lattice. If the band contains two electrons, 
the kinetic energy of each electron is —2t, the bottom of the band, hence 
Ekin = (In a tight-binding picture the reference energy is the center 
of the band irrespective of Ep, causing E^in to be always negative). Let 
us now consider the kinetic energy of a pair in the extreme pairing limit, 
i.e. J7 3> t, causing both electrons to occupy the same site, with an 
interaction energy —U. The occupation function in this case becomes 

This implies that the kinetic energy approaches E^in —St^ /U . Hence 
the kinetic energy increases from E'^^^ = —4t to Ef.-^ = — ^ when 
the local pairs are formed. The paired electrons behave like bosons of 
charge 2e. A second order perturbation calculation yields an effective 
boson hopping parameter[39] t' = t'^/U. In experiments probing the 
charge dynamics, this hopping parameter determines the inertia of the 
charges in an accelerating field. As a result the plasma frequency of such 
a model would be 

a;2 =4 "(26)2^^ (1.50) 
whereas if these pair correlations are muted 

aH 

^l,n = 47rne2^ (1.51) 

Because the plasma frequency is just the low frequency spectral weight 
associated with the charge carriers, this demonstrates, that for conven- 
tional pairs (i.e. those which are formed due to interaction energy lower- 
ing) the expected trend is, that in the superconducting state the spectral 
weight decreases. Note, that this argument can only demonstrate the 
direction in which the plasma frequency changes when the pair correla- 
tions become reduced, but it does not correctly provide the quantitative 
size of the change, since the strong coupling regime of Eq. 1.50 implies 
the presence of a finite fraction of uncondensed 'preformed' pairs in the 
normal state. The same effect exists in the limit of weak pairing corre- 
lations. In Ref. [40] (Eq. 29, ignoring particle-hole asymmetric terms) 
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the following expression was derived for the plasma resonance 

(1.52) 



dk 



where V is the volume of the system, and -E'l = + l^feP- Integrating 
in parts, using that A^E^^dkCk = dk (ek/Ek), and that dkCk = at the 
zone-boundary, we obtain 



a;.. = ^^2^— (1.53) 



f'^ V ^mu 

where mj.^ = h^'^d^€k/dk^ , and = 1 — ek/Ek- For a monotonous band 
dispersion the plasma frequency of the superconductor is always smaller 
than that of the unpaired system: Because the sign of the band-mass 
changes from positive near the bottom of the band to negative near the 
top, the effect of the broadened occupation factors is to give a slightly 
smaller average over m^^^, hence LOp is smaller. Note that the mass of 
free electrons does not depend on momentum, hence in free space LOp is 
unaffected by the pairing. 

To obtain an estimate of the order of magnitude of the change of 
spectral weight, we consider a square band of width W with a Fermi 
energy Ep = N^/ (21^), where Ne is the number of electrons per unit cell. 
To simplify matters we assume that 1 /m^ varies linearly as a function of 
band energy:l/m(e) = {W — lEp — 2e)/(l^mo). We consider the limit 
where A << W^Ep- Let us assume that the bandwidth ~ 1 eV, and 
A 14 meV corresponding to Tc =90 K. The reduction of the spectral 
weight is then 0.28 %. If we assume that the bandwidth is 0.1 eV, the 
spectral weight reduction would typically be 11.4 %. 



3.4 Kinetic energy driven superconductivity 

If the state above Tc is not a Fermi liquid, the situation could be 
reversed. The right-hand panel of Fig. 1.5 represents a state very dif- 
ferent from a Fermi liquid, and in fact looks similar to a gapped state. 
Indeed even for the ID Luttinger liquid n{k) has an infinite slope at 
kp- If indeed the normal state would have a broad momentum distri- 
bution like the one indicated, the total kinetic energy becomes lower 
once pairs are formed, provided that the slope of n{k) at kp is steeper 
in the superconducting state. This is not necessarily in contradiction 
with the virial theorem, even though ultimately all relevant interactions 
(including electron-phonon interactions) are derived from the Coulomb 
interaction: The superconducting correlations involve the low energy 
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scale quasi-particle excitations and their interactions. These effective 
interactions usually have characteristics quite different from the origi- 
nal Coulomb interaction, resulting in Ec/Ekin 7^—2 for the low energy 
quasi-particles. Various models have been recently proposed involving 
pairing due to a reduction of kinetic energy. In strongly anisotropic mate- 
rials such as the cuprates, two possible types of kinetic energy should be 
distinguished: Perpendicular to the plancs[41, 42] (along the c-direction) 
and along the planar directions [43, 44, 45, 46, 47, 48, 49, 50, 51]. 

4. Experimental studies of superconductivity 
induced spectral weight transfer 

4.1 Josephson plasmons and c-axis kinetic 

energy 

C-axis kinetic energy driven superconductivity has been proposed 
within the context of inter-layer tunneling, and has been extensively 
discussed in a large number of papers[43, 44, 41, 52, 42, 21, 53, 54, 55, 
56, 57, 20, 58, 59, 60, 61, 62, 63, 64]. One of the main reasons to suspect 
that superconductivity was c-axis kinetic driven, was the observation of 
"incoherent" c-axis transport of quasi-particles in the normal state [65] 
and, rather surprisingly, also in the superconducting statc[66, 67, 68], 
thus providing a channel for kinetic energy lowering for charge carriers 
as soon as pairing sets in. As discussed in section 2.4 a very useful tool 
in the discussion of kinetic energy is the low frequency spectral weight 
associated with the charge carriers. In infrared spectra this spectral 
weight is contained within a the 'Drude' conductivity peak centered at 
a; = 0. Within the context of the tight-binding model a simple relation 
exists between the kinetic energy per site, with volume per site T^, and 
the low frequency spectral weight [18, 19] 

Here the plasma frequency, uip, is used to quantify the low frequency 
spectral weight: 

uP' f^-rn 1 

+ / Rca{io)dio = -ujI (1.55) 
8 J0+ 8 

where the integration should be carried out over all transitions within 

the band, including the (5-function at a; = in the superconducting state. 

The 5{u)) peak in Recr(a;) is of course not visible in the spectra directly. 

However the presence of the superfluid is manifested prominently in the 

London term of Ree(a;) (proportional to lm(j(a;)): e/,(c<;) = — w^^^o;"^. 

In La2-j;Sr3;Cu04 the London term is manifested in a spectacular way 
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as a prominent plasma resonance perpendicular to the superconducting 
planes[69]. This is commonly used to determine the superfluid spec- 
tral weight, uJp g, from the experimental spectra. Apart from universal 
prefactors, the amount of spectral weight of the 6{lo) conductivity peak 
corresponds to the Josephson coupling energy, which in turn is the inter- 
layer pairhopping amplitude. It therefore provides an upper limit to the 
change of kinetic energy between the normal and superconducting state 
[41, 52], because the spectral weight transferred from higher frequencies 
to the 6{u)-peak cannot exceed this amount. This allowed a simple ex- 
perimental way to test the idea of c-axis kinetic energy driven supercon- 
ductivity by comparing the experimentally measured values of the con- 
densation energy {Econd) and Ej. The inter-layer tunneling hypothesis 
required, that Ej ~ Econd- In the spring of 1996 the first experimental 
results were presented[53] for T12201 (Tc=80 K), showing that Ej was 
at least two orders of magnitude too small to account for the condensa- 
tion energy. Later measurements of Ac [57] (approximately 17 jum) and 
the Josephson plasma resonance (JPR)[56] at 28 cm~^, allowed a defi- 
nite determination of the Josephson coupling energy of this compound, 
indicating that Ej ^ 0.3//eV in T12201 with = 80 K (see Fig. 1.7). 
This is a factor 400 lower than Econd ~ 100/^eV per copper, based either 
on cy experimental data[70], or on the formula Econd = 0.5A^(0)A^ with 
A^(0) = leV~^ per copper, and A ~ 15meV. In Fig. 1.8 the change in 
c-axis kinetic energy and the Josephson coupling energies are compared 
to the condensation energy for a large number of high Tc cupratcs. For 
most materials we see that Ej < Econd, sometimes differing by several 
orders of magnitude. 

These arguments falsifying the inter-layer tunneling mechanism have 
been questioned [21], arguing that a large part of the specific heat of 
T12201 is due to 3D fluctuations, and that these fluctuations should 
be subtracted when the condensation energy is calculated. However, 
it was recently shown [71] that due to thermodynamical constraints the 
fluctuation correction can not exceed a factor 2.5 in the case of T12201 
(as compared to a factor 40 in Ref. [21]). Hence the discrepancy between 
the Josephson coupling energy and the condensation energy of T12201 
is still two orders of magnitude. 

However, as stressed above, Ej provides only an upper limit loi AEj^i^- 
A c-axis kinetic energy change smaller than Ej is obtained if we take into 
account the fact that a substantial part of 5(a;)-function is just the spec- 
tral weight removed from the sub-gap region of the optical conductivity. 
Usually it is believed that in fact the latter is the only source of intensity 
of spectral weight for the (5-function, known as the (phenomenological) 
Glover- Tinkham-Ferell[72] sum rule. According to the arguments given 
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Figure 1.6. Grazing reflectivity of a Tl2Ba2Cu06 thin film (upper panel) and 
Lai.85Sro.i5Cu04 single crystal (lower panel) measured with the polarization of the 
incident light tilted at an angle of 80° relative to the copper-oxygen planes. For LSCO 
the Josephson plasma resonance can be clearly seen at 40 cm~^. For T12212 no the 
Josephson plasma resonance is observed, indicating that it is located below the lower 
limit of 30 cm~^ of the spectrometer. This implies that the Josephson coupling en- 
ergy in this compound is at least two orders of magnitude lower than required by the 
inter-layer tunneling hypothesis. Data from Ref. [54] 

in section 3.4 we may conclude that Ekin,.s = Ekin,n when we observe, 
that all spectral weight origins from the far-infrared gap region in agree- 
ment with the Glover-Tinkham-Ferrell sum rule. If, on the other hand, 
superconductivity is accompanied by a lowering of c-axis kinetic energy, 
part of uj'^^g originates from the higher frequency region of inter-band 
transitions, which begins at typically 2 eV. In other words, we may say 
that u^ g is an upper limit to the kinetic energy change 
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Figure 1.7. P-polarizcd reflectivity at 80° angle of incidence of Tl2Ba2Cu06. Prom 
top to bottom: 4K, 10 K, 20 K, 30 K, 40 K, 50 K, 60 K, 75 K, and 90 K. The curves 
have been given incremental 3 percent vertical offsets for clarity. Data from Ref. [56] 



A direct determination of Ekin,s — Ekin,n is obtained by measuring ex- 
perimentally the amount of spectral weight transferred to the 5{uj) peak 
due to the passage from the normal to the superconducting state, as 
was done by Basov et aZ.[20, 60]. These data indicated that for under- 
doped materials about 60% comes from the sub-gap region in the far 
infrared, while about 40% originates from frequencies much higher than 
the gap, whereas for optimally doped cuprates at least 90% originates 
from the gap-region, while less than 10% comes from higher energy. Ex- 
perimental artifacts caused by a very small amount of mixing of ab-plane 
reflectivity into the c-axis reflectivity curves may have resulted in an 
overestimation of the spectral weight originating from high energies [60], 
in particular those samples where the electronic cJc(^) is very low due 
to the 2-dimensionality. Optimally doped YBCO is probably less prone 
to systematic errors due to leakage of Rah into the c-axis reflectivity, 
since aduj) of this material is among the largest in the cuprate family. 
The larger (Tc(w) causes the c-axis reflectivity to be much larger at all 
frequencies, thereby reducing the effect of spurious mixing of a6-plane 
reflectivity in the optical spectra on the Kramers-Kronig analyzes. 

In summary AEkin,c < O.lEj in most cases. For several of the single- 
layer cuprates it has become clear now, that AEkin significantly under- 
shoots the condensation energy, sometimes by two orders of magnitude 
or worse, as indicated in Fig. 1.8. 
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Figure 1.8. Intrinsic Josephson coupling energy [20, 55, 56, 57, 58, 59, 66, 69, 75] 
versus condensation energy[70, 76] 



4.2 Josephson plasmons in mult i- layered 
cuprates 

This situation may be different for tfie bi-layer compounds. In these 
materials in principle the coupling within the bi-layer may provide an 
additional source of frustrated inter-layer kinetic energy, which can in 
principle be released when the material enters the superconductng state. 
This can in principle be monitored with infrared spectroscopy, because 
quite generally a stack of Josephson coupled layers with two different 
types of weak links alternating (in the present context corresponding to 
inter-bilayer and intra-bilayer) should exhibit three Josephson collective 
modes instead of one: Two of those modes are longitudinal Josephson 
plasma resonances, which show up as peaks in the energy loss function 
\Ta.{—\/e{u))). In between these two longitudinal resonances one expects 
a transverse optical plasma resonance, which is revealed by a peak in 
Recj(ci;). In essence the extra two modes are out-of-phase oscillations 
of the two types of junctions. This has been predicted in Ref. [77] for 
the case of a multi-layer of Josephson coupled 2D superconducting lay- 
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ers. Further detailed calculations for the bi-layer case were presented 
in Refs. [78, 79]. The existence of two longitudinal modes and one 
associated transverse plasmon mode at finite frequencies has been con- 
firmed experimentally for the SmLao.8Sro.2Cu04_5 in a series of papers 
[80, 81, 82, 83, 84] (see Fig. 1.9). 

The c-axis optical conductivity of YBCO is one order of magnitude 
larger than for LSCO near optimal doping. As a result the relative im- 
portance of the optical phonons in the spectra is diminished. In the case 
of optimally doped YBCO, the experiments indicate no significant trans- 
fer of spectral weight from high frequencies associated with the onset of 
superconductivity. C-axis reflectivity data[75] of optimally doped YBCO 
are shown in Figs. 1.10. Above Tc the optical conductivity is weakly fre- 
quency dependent, and docs not resemble a Drude peak. Below Tc the 
conductivity is depleted for frequencies below 500 cm~^, reminiscent of 
the opening of a large gap, but not an s-wave gap, since a relatively large 
conductivity remains in this range. 

There is a slight overshoot of the optical conductivity in the region 
between 500 and 700 cm~^, due to the fact that the normal state and 
superconducting state curves cross at 600 cm~^. In the case of the bi- 
layer cuprates this could be explained as a result of the presence of two 
superconducting layers per unit cell, resulting in the 'transverse optical' 
plasma mode mentioned above[62, 63, 64, 73, 74, 75]. 

For the f-sum rule the presence of this extra mode makes no difference. 
The extra spectral weight in the superconducting state associated with 
this mode has in principle the same origin as the spectral weight in the 
zero-frequency 5-function. In a conventional picture the source would 
be the spectral weight, removed due to a depletion of adu}) in the gap- 
region. The implementation of the sum rule relevant for this case then 
states that the relative spectral weight function 

overshoots the 100 % line close to the 'second plasma' mode, and sat- 
urates at 100 % for frequencies far above this mode. This is indeed 
observed in Fig. 1.10. 

Additional studies of the bi-layer (and tri-layer) materials have pro- 
vided conflrmation of the transverse optical plasmon in these materials. 
In spite of its high frequency, making the assignment to the Josephson 
effect rather dubious, nevertheless the transverse optical mode either 
makes its first appearance below Tc, or gains in sharpness and intensity 
at the temperature where pairs are being formed (which for under-doped 
cuprates begins already above Tc). Also in at least a number of cases 
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Figure 1.9. (a) Real part of the c-a:xis dielectric function of SmLao.8Sro.2Cu04_5 for 
4 K (closed symbols), and 20 K (open symbols) (b) The c-axis loss function, Ime(a;)~^. 
(c) Real part of the c-axis optical conductivity. Data from Ref. [83, 85]. 



the spectral weight of the 'transverse optical' plasmon observed below 
Tc appears to originate not from the spectral weight removed from the 
gap region, but from much higher energies [60, 61, 62, 63, 64]. The impli- 
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Figure 1.10. C-axis optical spectra of optimally doped (x=6.93) and over-doped 

(x=7.0) YBa2Cu307-a;. From top to bottom: reflectivity, optical conductivity, 
dynamical impedance and relative spectral weight (Eq. 1.57). The dynamical 
impedance, pi(a;)=Re47r/ci;e(aj) is proportional to the energy loss function weighted 
by a factor l/cj. The optical phonons have been subtracted from the loss- functions 
for clarity. The data are from Ref. [75, 86] 



cation of this may be, that a non-neghgible fraction of frustrated c-axis 
kinetic energy is released when these materials become superconduct- 
ing. This seems to be particularly relevant for the strong intra-bilayer 
(or tri-layer) coupling of Bi2212, Bi2223 and Y123. 
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4.3 



Kinetic energy parallel to the planes 



In-plane kinetic energy driven superconductivity has been proposed 
by a number of researchers: Hirsch[43, 44, 87] discussed this possibiHty 
as a consequence of particle-hole asymmetry. It has also been discussed 
within the context of holes moving in an anti-ferromagnetic background 
[88, 89, 90, 91, 92, 93]. More recently the possibility of a reduction of 
kinetic energy associated with pair hopping between stripes has been 
suggcstcd[47, 51], and an in-plane pair dclocalization mechanism have 
been proposed in the context of the resonating valence band model[49. 



A major issue is the question how to measure this. The logical ap- 
proach would be to measure again a{uj,T) using the combination of 
reflectivity and Kramers-Kronig analysis, and then compare the spec- 
tral weight function in the superconducting state to the same above Tc. 
There are several weak points to this type of analysis. In the first place 
there is the problem of sensitivity and progression of experimental er- 
rors: Let us assume, that the change of kinetic energy is of order 0.1 
meV per Cu atom (this is approximately the condensation energy of the 
optimally doped single layer cuprate T12201, with Tc = 85 K.). For 
an inter-layer spacing of 1.2 nm, this corresponds to a spectral weight 
change A(t'p) = 10^ cm~^. As the total spectral weight in the far in- 
frared range is of order ly^ = 14000^ cm , the relative change in spectral 
weight is of order 0.05 %. Typical accuracy reached for spectral weight 
estimates using conventional reflection techniques is of order 5%. This 
illustrates the technical difficulties one has to face when attempting to 
extract superconductivity induced changes of the kinetic energy. 

Experimental limitations on the accuracy are imposed by (i) the im- 
possibility to measure all frequencies including the sub-mm range, (ii) 
systematic errors induced by Kramers-Kronig analysis: The usual pro- 
cedure is to use data into the visible/ultra-violet range and beyond for 
completing the Kramers-Kronig analysis in the far infrared, assuming 
that no important temperature dependence is present outside the far 
infrared range. Obviously this assumption becomes highly suspicious if 
the search is concentrated on spectral weight transfer originating from 
precisely this frequency range. 

The remedy is, to let nature perform the spectral weight integral. Due 
to causality Ree(a;) and Rec7(a;) satisfy the Kramers-Kronig relation 



The main idea of spectral weight transfer is, that spectral weight is 
essentially transferred from the inter-band transitions at an energy of 



50]. 




(1.58) 
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several eV, down to the 5-function in <t(u;) at u; = 0. Indeed various 
groups have reported a change of optical properties in the visible part of 
the spectrum when the sample becomes superconducting[96, 97, 98, 99]. 
If this is the case, we have a; = for the extra spectral weight in relation 
1.58. Together with Eq. 1.56 it follows that changes in kinetic energy 
can be read directly from Ree(a;) using the relation 



^E^li^) = -^^e6e{u) (1.59) 



If the spectral weight is transferred to a frequency range uiq, then the 
above expression can still be applied for a; » o^q. If we measure Ree(a;) 
directly using spectroscopic ellipsometry, then indeed nature does the 
integration of (y{ijo) for us at each temperature. This eliminates to a 
large extent various systematic errors affecting the overall accuracy of the 
spectral weight sum. It is important to measure the complex dielectric 
constant for a large range of different frequencies. 

The second problem is that already above the superconducting phase 
transition the optical spectra of these materials have appreciable tem- 
perature dependence. What we really like to measure is the spectra 
of the same material in the superconducting state, and in the 'normal' 
state, both at the same temperature. Typical magnetic fields required 
to bring the material in the normal state are impractical, let alone the 
complications of magneto-optics which then have to be faced. A more 
practical approach is to measure carefully the temperature dependence 
over a large temperature range, with small temperature intercepts, and 
to search for changes which occur at the phase transition. 

In Fig. 1.11 the spectral weight from to 10000 cm~^ is shown as a 
function of temperature for the case of Bi2212[99]. Note that this inte- 
gral corresponds to minus the ab-plane kinetic energy. We observe, that 
in the superconducting state the kinetic energy drops by an amount of 
about 1 meV per Cu. This is in fact a relatively large effect. This surpris- 
ing result seems to tell us that in the cuprates the kinetic energy in the 
superconducting state is lowered relative to the normal state. This corre- 
sponds to the unconventional scenario depicted in the right-hand panel 
of Fig.1.5, where the normal state is a non Fermi liquid, whereas the su- 
perconducting state follows the behavior of a (more) conventional BCS 
type wave-function with the usual type of Bogoliubov quasi-particles. 
The amazing conclusion from this would be, that there is no need for a 
lowering of the interaction energy any more. The condensation energy 
of optimally doped Bi2212 is about 0.1 meV per Cu atom[76]. 



Optical signatures of electroncorrelations in the cuprates 



33 



<N 
> 



Q 

+ 



< 4.10 



00 



> 

+ 



< 3.60 



CM 

00 




(10^ K^) 



Figure 1.11. Measured values of the quantity c^ujp^s + 8 Re(Tai,(a»)dai of Bi2212 
(Tc=88 K). The data are taken from Ref. [99, 100]. To make the conversion to kinetic 

energy summed over the two afo-planc directions, the numbers along the vertical axis 
have to be multiplied with a factor -10^K/(47re^a^) = -83 meV / eV^. 



5. Conclusions 

The optical conductivity is a fundamental property of solids, contains 
contributions of vibrational and electronic character. Among the elec- 
tronic type of excitations the intra-band and inter-band transitions, exci- 
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tons, and plasmons of different types correspond to the most prominent 
features in the spectra. In addition multi-magnon excitations or more 
exotic collective modes can often be detected. The careful study of the 
optical properties of solids can provide valuable microscopic informa- 
tion about the electronic structure of solids. In contrast to many other 
spectroscopic techniques, it is relative easy to obtain reliable absolute 
values of the optical conductivity. As a result sum rules and sum rule 
related integral expressions can often be applied to the optical spectra. 
Here we have treated a few examples of sum rule analysis: Applica- 
tion of the f-sum rule to the phonon spectra of transition metal silicides 
provides information on the resonant electron-phonon coupling in these 
materials. Integration of the energy-loss function gives the value of the 
Coulomb energy stored in the material, which is seen to increase when a 
high Tc cupratc enters the superconducting state. The spectral weight 
within the partially band of the high Tg cuprates is seen to become larger 
in superconducting state. This effect exists both perpendicular to the 
planes and parallel to the planes. This spectral weight change can be 
associated with a decrease of kinetic energy when the material becomes 
superconducting. Although the relative spectral weight change along 
the ab-plane is quite small, it indicates a fairly large change of the ab- 
plane kinetic energy, large enough to account for the energy by which the 
superconducting state of these materials is stabilized. In addition the 
real and imaginary part of the optical conductivity can be used to study 
the intrinsic Josephson coupling between the superconducting planes. 
In superconductors with two or more different types of weak links al- 
ternating, such as SmLai_a;Sra;Cu4, YBa2Cu307_a;, Bi2Sr2CaCu208, a 
rich spectrum of plasma-oscillations is observed in the superconducting 
state, and sometimes above Tc, due to the multi-layered structure of 
these materials. This has provided important insights in the nature of 
the coupling, and it has been used to extract quantitative values of this 
coupling. 
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